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Abstract 

The correspondence between classical and quantum invariants is established. The Ermakov Lewis 
quantum invariant of the time dependent harmonic oscillator is translated from the coordinate and 
momentum operators into amplitude and phase operators. In doing so, Turski's phase operator 
as well as Susskind-Glogower operators are generalized to the time dependent harmonic oscillator 
case. A quantum derivation of the Manley-Rowe relations is shown as an example. 
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I. INTRODUCTION 



Exact invariants have been extensively used to solve the time dependent Shrodinger 
equation [1]. Various related invariants have been obtained for the quantum mechanical time 
dependent harmonic oscillator equation in one dimension (QM-TDHO) . The Ermakov Lewis 
invariant and orthogonal functions invariant are two such constants of motion that have been 
used to solve the QM-TDHO problem. The Ermakov Lewis invariant is usually expressed in 
terms of coordinate and momentum operators although it has also been expressed in terms 
of raising and lowering operators that lead to number states for wave functions that are 
eigenstates of the invariant operator [2]. However, the amplitude operator stemming from 
this procedure does not correspond to the amplitude of the oscillator in the classical limit. 

The purpose of this communication is to translate the invariant formalism from the 
coordinate and momentum operators into an invariant in terms of amplitude and phase 
operators that reduce to the corresponding variables in the classical limit. In the second 
section, the solution to the QM-TDHO equation is stated using the square of the orthogonal 
functions invariant and the Ermakov invariant. In the following section, a second linear 
Hermitian invariant is introduced and the Ermakov-Lewis invariant is economically obtained 
from these two constants of motion. Two distinct annihilation and creation operators are 
presented in section four and their equations of motion are established. In section five, the 
quantum phase is defined using the Turski and Susskind and Glogower formalisms. The 
former definition is shown to yield an amplitude and phase representation that is consistent 
with the classical limit. In the last section, the Ermakov Lewis invariant is written in 
amplitude and phase variables. The energy conservation and photon number relations in 
nonlinear optical processes is shown as an example. 

II. EVOLUTION OPERATORS AND INVARIANTS 

Consider the time dependent Schrodinger equation with h = 1 

>^ ^ mm). (1) 

The solution to this equation for a time-independent Hamiltonian is formally given by 
\ip{t)) = U{t)\il){0)) , where U{t) is the evolution operator U{t) = exp l—iHt). For the time 
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dependent harmonic oscillator Hamiltonian 

m = \{p'+^\t)cf). (2) 

the solution may be written in terms of a propagator that involves a time independent 
operator together with an appropriate transformation of the wave function 

m))-uiT^m\m)- (3) 

The propagator is given by 

"* dt' 



Uj = exp (^-iSaL^ ; Sa = J (4) 



where the function Sa is a time dependent c-number and the transformation is defined as 

.In a dq'^\ ^ f .d\na(p\ ^ f .\\ia^^^ ^ _ ( . ol ^2 



T = exp j cxp ) = exp [i^{qp + pq) j exp j . (5) 

The time independent operator in the propagator is an invariant that is not unique [3]. 
On the one hand, it may be proportional to the square of the orthogonal functions invariant 
operator 

^« = ^ {up - uqf , (6) 

where the function a — > e M, replaced in the invariant as well as in the transformation 
expressions, obeys the TDHO equation 

u + n'^{t)u = Q. (7) 

On the other hand, the propagator may be written using the Ermakov Lewis invariant 
with a — > p, where p obeys the Ermakov equation 

p + Q'{t)p^p-\ (8) 

In either case, it is seen that the invariant in the time dependent case enters the propagator 
expression in an analogous fashion as the Hamiltonian does in the time independent case. 

III. CLASSICAL AND QUANTUM INVARIANTS 

The classical orthogonal functions invariant is 

G = ?ig2-?2?i, (9) 
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where qi and q2 are real linearly independent solutions of the TDHO equation [4]. The 
quantum invariant arising from the mapping of q2 and q2 into the coordinate and momentum 
operators is 

Gi = Uip — iiiq. (10) 

The obtention of a second invariant given a first invariant has been subject of several com- 
munications [5] , [6] . It is worth remarking that the existence of a second invariant warrants 
complete integrability for a Hamiltonian Ermakov system [7] . Within the present formalism, 
it is straight forward to introduce a second invariant stemming from the mapping of qi and 
qi into the coordinate and momentum operators 

G2 = -U2P + U2q- (11) 

These two invariants obey the commutation relation [Gi, G2] = —iG, where G is a con- 
stant. From the sum of their squares, we may construct the invariant operator 

I-l{Gl + Gl), (12) 
which in terms of the position and momentum operators is 

Gl + Gj^ {uip - iiiqf -\- (-U2P + U2qf 

= {u\ + u\)p^ + {uf + ul)(f - {uiiii + U2U2){pq + qp)- (13) 



We may rewrite this expression as a function of an amplitude function p = \J\i\ + by 
noticing that pp = Uiiii + U2U2 and that the orthogonal functions obey (9), so that 

.2 {UiU2 - U2Ulf + {Uiili + U2U2f / ■ 2 , ■ 2\ M a\ 

-J + p = = (Ml + M2) • 14 

The operator in terms of p is then 

(15) 



7^ 
2 



^ ) + (pp - 



but this is precisely the Ermakov Lewis invariant where the real constant G is usually 
normalized to unity in general different from one [8]. The above procedure is a simple 
derivation of the quantum Ermakov Lewis invariant, which has otherwise been obtained 
using rather more complex mathematical methods [9]. The non Hermitian linear invariant 
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7c introduced by Malkin et al. written in terms of the orthogonal functions invariants is 
7c = Gi - iG2- 

In the classical case, the amplitude and phase representation of the invariant is straight 
forward. If the coordinate variable u is expressed in polar coordinates 

pe''" ^ape-''", (16) 

where cr is a constant, the classical orthogonal functions invariant (9) may be expressed in 
terms of amplitude and phase variables as 

Gl (1 - a') = /i,. (17) 

The constant G/ (1 — cr^) may be normahzed to one and the derivative of the phase 
written as the frequency uj {t) = s^; the squared amplitude times the frequency then obey 
the relationship 

p^uj{t) = l. (18) 

The energy of a time independent oscillator is proportional to the squares of the mo- 
mentum and coordinate variables £ oc -\- cu^q^, which in terms of the amplitude and 
phase variables is E (x PqOJq. If this relationship is considered to hold in the time depen- 
dent case, the invariant is then proportional to the ratio of the energy over the frequency 
G/ (1 — a^) cx £ {t) /u) {t) thus yielding the well known adiabatic invariant [10]. Therefore, 
the classical orthogonal functions invariant is proportional to the eigenvalue of the Hamil- 
tonian function of the system. 

Nonetheless, the quantum versions of this invariant produces a linear form in the coor- 
dinate and momentum operators as seen in Eqs. (10) and (11). It therefore comes to no 
surprise that the argument of the propagator is proportional to the square of the orthogonal 
functions quantum invariant (6). On the other hand, the classical Ermakov Lewis invari- 
ant in the amplitude and phase representation follows from the substitution q — > p cos Sp, 
p^dq/dt [8]: 

I = \p'sl (19) 

that implies a quadratic dependence on the energy of the oscillator. Therefore, a quantum 
invariant with a quadratic dependence on the coordinate and momentum variables should 
be in correspondence with the classical orthogonal functions invariant. 
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IV. CREATION AND ANNIHILATION OPERATORS 



An operator that can be written as the sum of two squares may be expressed in terms of 
two adjoint complex quantities. To wit, given an operator (5 that can be expressed as 

/9 = fc? + fe2, (20) 

provided [&i,&2] = c, with c a c-number, there exist annihilation and creation operators 
h — hi-\- ib2, — bl — ib2 so that the operator may be written as P — b^b — i 61, 62 ■ For 
instance, annihilation and creation operators for the Hamiltonian (2) may be written as [11] 

{n'/'{t)q + tPin''\t)) , = ^ {^''\t)q - tP/n''\t)) . (21) 

However, the way in which the $ operator is written as the sum of two squares need not 
be unique. In fact, for the invariant operator defined in the previous section, expressions (12) 
and (15) are two such possibilities. The former leads to annihilation and creation operators 
of the form 

i = ^ (Gi - tG2) , i1" = (Gi + 1G2) , (22) 

where the identification 61 Gi and 62 —G2 has been made. These operators may also 
be obtained from the non Hermitian linear invariant which arise from the complex solution 
of the TDHO equation [12], [13]. The sign in the imaginary part of the above expressions 
is introduced in order to have consistency with the cited results. These annihilation and 
creation operators are also invariant since they are composed by invariant operators. 
On the other hand, the operators arising from (15) yield 

« (^) = ^ + i{PP - PQ?) , «^ it) = 7f - '(PP - P^^) ■ (23) 

These time dependent annihilation and creation operators were originally introduced by 
Lewis [2] . The Ermakov invariant in terms of these operators is 

i^dHt)d{t) + l^A^A + ^. (24) 

In order to obtain the transformation between the distinct annihilation and creation opera- 
tors, evaluate 

^is,i^^-is,i ^ ^^-is, ^A_(g^- iG^^ (^U2 + iui) (25) 
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where ui — —psinsp, U2 — p cos Sp. The relationships between the orthogonal functions 
and their trigonometric representation is not unique. This choice represents the function ui 
leading U2 by | as Sp increases [10]. Replacing the definitions of the invariants yields 
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U1U2 - UiU2\ . , . /m? + . / uiiii + U2U2 , „ 
~P )'^\~^)''~\ ~P 



(26) 



which simplifies to 

ie-^ = -L + ^{pp - pq)^ = a. (27) 

Therefore the time dependent annihilation (creation) operators may be written as the prod- 
uct of the time independent annihilation (creation) operators times a phase that only involves 
a c-number function. This expression may be written as a unitary transformation of a phase 
shift 

a — exp (^ispi^ A exp ^—ispl^ (28) 
The equation of motion of this operator is then 

h^iu;{t)[i,a\. (29) 

It is thus seen that the operator uj{t)I in the QM-TDHO again plays the role that the 
Hamiltonian does in a time independent harmonic oscillator case. This assertion is consistent 
with the transformation that relates the invariant and the time dependent Hamiltonian [3] 

u;{t)I = Hit) - z-^T. (30) 

V. PHASE OPERATORS FOR TIME DEPENDENT OSCILLATOR 

As it is well known, different operators can be used to define the phase in quantum optics 
[14]. The invariant formalism will be applied here to the phase operator given by Turski 
and the exponential phase operators of Susskind and Glogower. In particular, the former 
operator will allow an appropriate translation of the classical amplitude-phase invariant into 
the quantum one. 



A. Turski phase operator 

By using the annihilation operator (23) the displacement operator can be written as 
D{a) — exp {aa^ — a*a), a — rexjp{i9). The vacuum state may then be displaced to obtain 
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a coherent state \a) — £>(q;)|0) and the phase operator introduced by Turski [15] is then 
generahzed to the time dependent case 



$ = y e\a){a\d^a. (31) 

This operator obeys the commutation relation /] = —i. In order to evaluate the time 
evolution of this operator can be written in terms of the invariant annihilation and 

creation operators using (28) 

$ = e''"^ (^J eDA{a)e-''''^\0){0\e''''^ID\{a)(fa^ 6-''"^ (32) 

where Da{.(^) = exp {aA'^ — a* A). The invariant acting over the vacuum state is /|0) = ||0) 
and the phase is then 

$ = e''"' Q eDAia)\0){0\b\{a)d^a^ e-''"' , (33) 

the time derivative of this expression yields the equation of motion for 

^ ^iu;(t)[I,^^ -u;{t). (34) 

The operator u!{t)i once again takes the role of the Hamiltonian since the Turski operator 
commutes with H in the time independent case. 

B. Susskind-Glogower operators 

The generalization of the phase to the time dependent case is also applicable using other 
formalisms. Consider, for example, the Susskind-Glogower operators [16] given by (see for 
instance [17]) 



V 



^ oo ^ oo 

_a = V|n)(n + l|, 1/^ = a^^== = V |n + l)(n|, (35) 



where \n) is a number state, eigenstate of the invariant /. The unitary transformation 
VIV^ — /-|-l,works as a shifter in the same way as q andp do: exp{iap)q exp{—iap) — q+a. 
The sine and cosine operators for the Susskind-Glogower operators 



give the commutation relations [/, C] — —iS, [/, S] — iC. Following the same treatment as 
above, i.e. expressing operators that depend on a and in terms of the invariants A, 
and /, the equations of motion for the sine and cosine operators are 

C ^iuj{t)[i,C\^uj{t)S, S^iuj{t)[i,S\^-uj{t)C. (37) 



VI. AMPLITUDE AND PHASE REPRESENTATION OF INVARIANT 

The coordinate operator from (23) is 

and following Dirac [18] the creation and annihilation operators may be written as 

a = \/7e-'*, a) = e'^^Tl. (39) 
The coordinate operator (38) in the form of amplitude and phase variables is then 



2c; (t) \ 2uj{ty 

where the amplitude p and phase Sp are identified as 



^ e-^* + e**W-^, (40) 



The orthogonal functions invariant with the aid of (34) is given in amplitude and phase 
operators as 

which has the same structure of the orthogonal functions classical invariant written in am- 
plitude and phase variables (17). The number operator is then identified with 

^ = 43 

This identification may seem dimensionally awkward but it should be remembered that 
the invariant initial value was normalized to one (18). The explicit introduction of the nor- 
malization factor PqCUq makes of course a dimensionless photon number for an adimensional 
amplitude po- Since d^d is invariant from (24), if the frequency is constant the number of 
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excitations is then also constant. However, in the time dependent case, the number of ex- 
citations is inversely proportional to the time dependent frequency in correspondence with 
the intensity dependence obtained in the classical limit. 

The energy of the excitation at a given time tg is given hj S = nu {tg) = n<l> (with h = 1) 
and thus the quantum invariant represents the energy conservation of the system. Consider, 
as an example of this formalism, the number of excitations to represent the photon number. 
Allow for a non degenerate nonlinear process where the field experiences second harmonic 
generation (SHG). Let the photon number at time ti be rii when the frequency mode is cui 
and allow it to evolve at a time t2 to the mode UJ2 — 2a;i. The number of photons in the 
mode U2 is then 

712 = ni — = -ni. (44) 

This scheme corresponds to a Lagrangian hydrodynamic framework where a given volume 
is being followed along its propagating path. On the other hand, an Eulerian framework 
in a fixed point in space under steady state conditions corresponds to the cigenfunctions 
of the invariant operator (24). The photon number is then a^a and is constant for each 
participating mode. If the invariant is the same for all modes, the energy of the modes are 
related by 

- = -• (45) 

In order to preserve the total energy of the system two degenerate modes with energy 
£1 are of course required in the SHG case. This reasoning may be extended to an arbitrary 
number of modes leading to other nonlinear processes such as parametric amplification or 
frequency difference. These type of equations are known in nonlinear optics as Manley-Rowe 
relations [19]. They are usually derived in semiclassical theory through a rather cumber- 
some procedure that relies on the particular nonlinearity being described together with the 
symmetries that they involve (Kleinmann's condition) when no absorption is present [20]. 
These semiclassical results are often interpreted in terms of photon numbers participating 
in each mode [20], [21]. This interpretation is naturally embodied in the present quantum 
treatment. 
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VII. CONCLUSIONS 



An economical derivation of the quantum Ermakov Lewis invariant has been presented. 
This invariant may be used in an equivalent fashion as the Hamiltonian is used in the time 
independent case. Namely, to obtain evolution operators, to cast the equations of motion 
of different operators in commutative expressions, and to produce a phase shift with its 
exponential form. The invariant and time dependent definitions for annihilation and creation 
operators have been used to generalize the quantum phase to the time dependent case. 
Following the classical form of the orthogonal functions invariant, the quantum Ermakov 
Lewis invariant has been expressed in amplitude and phase variables in accordance with 
the correspondence principle. A quantum derivation of the Manley-Rowe relations has been 
presented as a particular application of this representation. 
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